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ABSTRACT 

Following re-analysis of Hubble Space Telescope observations of primary transits of the extrasolar 
planet HD 209458b at Lyman a, Ben-Jaffel (2007, BJ007) claims that no sign of evaporation is ob- 
served. Here we show that, in fact, this new analysis is consistent with the one of Vidal-Madjar et 
al. (2003, VM003) and supports the detection of evaporation. The apparent disagreement is mainly 
due to the disparate wavelength ranges that are used to derive the transit absorption depth. VM003 
derives a (15 ± 4)% absorption depth during transit over the core of the stellar Lyman a line (from 
— 130 km/s to +100 km/s), and this result agrees with the (8.9 ± 2.1)% absorption depth reported 
by BJ007 from a slightly expanded dataset but over a larger wavelength range (±200 km/s). These 
measurements agree also with the (5 ±2)% absorption reported by Vidal-Madjar et al. (2004) over the 
whole Lyman a line from independent, lower-resolution data. We show that stellar Lyman a variabil- 
ity is unlikely to significantly affect those detections. The Hi atoms must necessarily have velocities 
above the escape velocities and/or be outside the Roche lobe, given the lobe shape and orientation. 
Absorption by H I in HD 209458b's atmosphere has thus been detected with different datasets, and 
now with independent analyses. All these results strengthen the concept of evaporating hot-Jupiters, 
as well as the modelization of this phenomenon. 

Subject headings: planetary systems — stars: individual (HD 209458) — spectroscopy 
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1. INTRODUCTION 

Few detections of extrasolar planets' atmospheric 
species are reported so far, but they have been recog- 
nized as important steps in our understanding of these 
objects. Of particular interest are the direct detections 
with the Hubble Space Telescope (HST) during primary 
transits of Nal by Charbonneau et al. (2002) as well as 
HI, OI and CII by Vidal-Madjar et al. (2003, 2004) (here- 
after VM003 and VM004) and Ballester et al. (2007). 

The recent paper by Ben Jaffcl (2007) (hereafter 
BJ007) however casts some doubt on many aspects of 
the HI detection in the upper atmosphere of HD 209458b 
and on the implication that the planet is evaporating 
due to the large energy input from its nearby host star 
(e.g., Lecavelier des Etangs 2007). Consequently, a re- 
lated large number of theoretical studies would have to 
be revised. 

In the present rebuttal paper, we discuss the BJ007 
arguments and show where they are misguiding. 

2. WHERE IS THE DIFFERENCE ? 

BJ007 completed a new data analysis based on sam- 
pling the same observations as VM003 in a different tem- 
poral manner, and adding data from two HST orbits 
from an archival program completed at another epoch 
(two orbits added to the nine HST orbits of the VM003 
observations). From the resulting new Lyman a tran- 
sit light curve, the transit absorption depth is evaluated 
by integrating the Lyman a flux in two ranges : "blue" 
(1214.83-1215.36A) and "red" (1215.89-1216.43A). This 

* Because the escape of hydrogen atoms is strengthened in this 
paper, we renew our proposal to use the nickname "Osiris" for the 
planet HD 209458b who loses mass like the Egyptian god. 



slightly increases the wavelength domain excluded from 
the analysis because of geocoronal contamination, from 
1215.5-1215.8A in VM003 to 1215.36-1215.89A in BJ007. 
This has, however, no significant consequences on the 
evaluation as shown in the Figure 4 of VM003. 

BJ007 evaluates the transit depth by fitting a light 
curve over the observations sampled by 300s sub- 
exposures as a function of the planet orbital phase. The 
resampled observations show some level of variability 
mainly due to the variation of the stellar Lyman a flux. 
The fitting procedure "smooths" these variations as if 
all observations were made with an average stellar Ly- 
man a flux. This gives an average HI absorption during 
the transit of (8.9 ± 2.1)% by considering a wavelength 
range from 1214.83 to 1216.43 A which corresponds to 
~ ±200 km/s in velocity space. 

The data analysis of BJ007 is not put into question 
here. The major differences between BJ007 and VM003 
are on the data interpretation. Both analyses provide a 
Lyman a line flux as a function of time and wavelength, 
which includes possible stellar variations and planetary 
transit signature. The differences lie i) in the wavelength 
ranges used to convert the spectra, as a function of time, 
into a single absorption depth measurement and ii) in 
the reference flux used to correct for the intrinsic stellar 
flux variations. 

The Lyman a emission of the star can significantly 
vary from epoch to epoch. Because there is no detectable 
transit absorption signature in the Lyman a line wings 
(the nominal 1.5% obscuration by the planetary disk is 
below the data S/N), VM003 calculated a relative ab- 
sorption depth using wings of the Lyman a line as flux 
reference. This method aims at correcting for any intrin- 
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Fig. 1. — The observed HD 209458 Lyman a profiles as observed 
by Vidal-Madjar et al. (2003) before and during the planetary tran- 
sit. The BJ007 re-analysis of nearly the same data set produced a 
similar Lyman a line profile. The two vertical dashed lines define 
the limits Ai and A2 of the line core where HI planetary absorp- 
tion takes place. In VM003 as well as in BJ007, the central part of 
the line (noted "Geo") possibly perturbed by the Earth geocoro- 
nal emission is omitted from the analysis. The line wings are used 
by VM003 as a flux reference to correct for the stellar Lyman a 
intrinsic variations. 

sic and unknown changes in the Lyman a stellar flux (see 
§ 13). As illustrated in Fig. [TJ VM003 defined two spec- 
tral domains : the line core from Ai to A2 (called "In" 
in VM003 and "Core" in Fig. [TJ excluding the part of 
the spectrum contaminated by the geocoronal Lyman a 
emission), and the remaining wavelength domain in the 
wings of the line used as a flux reference from I2I4.4A to 
Ai and from A 2 to 1216.8A (called "Out" in VM003 and 
"Wing" in Fig.QJ. The best domain defined by Ai and A2 
is found by varying the positions of these two wavelengths 
until the strongest absorption signal is identified. VM003 
found that the strongest absorption relative to the line 
wings takes place in the line core between Ai=1215.15 
and A 2 = 1216.10A (excluding the central geocoronal re- 
gion). 

Because the HI absorption takes place only in the cen- 
tral part of the Lyman a line, the absorption depth mea- 
surement must decrease for increasing wavelength range, 
including wings of the line where there is no absorption. 
This has already been found by VM004 with an inde- 
pendent data set obtained with STIS in the G140L low 
resolution mode. In the low resolution data, the Lyman a 
stellar emission line is not spectrally resolved and only 
the total Lyman a flux can be evaluated. VM004 ob- 
tained a (5 ±2)% transit absorption depth over the whole 
line, in agreement with the estimate obtained from the 
VM003 dataset (Table [J). 

To compare the VM003 result to that of BJ007, we 
can also calculate the transit absorption depth over the 
spectral domain as defined by BJ007 using the spectra 
of VM003. Keeping the same approach as in VM003 
to account for stellar Lyman a variations, we evaluate 
the planetary absorption during transit using the line 
wings as flux reference. We find a mid-transit absorp- 



tion of (7.3±2.0)% over the same spectral region used by 
BJ007, which is in agreement with the BJ007 result of 
8.9% considering that the data set is the same, except 
for the addition of two HST orbits to the nine used by 
VM003. This result shows that, considering the same 
wavelength range, similar transit absorption depths are 
found in the VM003 and BJ007 spectra (Table [T]). 

The uncertainty and noise appear to be lower in BJ007 
than in VM003. The difference is explained by the larger 
wavelength range used by BJ007 to estimate the Ly- 
man a flux and by the uncertainty introduced by the 
correction of the Lyman a variations applied by VM003 
(see§El). 

In short, the results given in VM003 and BJ007 for Ly- 
man a absorption depths can be reconciled. The appar- 
ent difference is basically due to the width of the spectral 
domain over which the absorption is computed, acknowl- 
edging that the HI absorption does not cover the whole 
extent of the stellar Lyman a line. When using a larger 
wavelength domain, the absorption signal is diluted and 
the absorption depth measurement is lower, as found in 
BJ007 compared to VM003, and in VM004 compared to 
BJ007. 

3. HIGH VELOCITY BLUE-SHIFTED ABSORPTIONS 

Another argument, made by BJ007 against the evapo- 
ration scenario, is that the blue-shifted absorption (pro- 
duced by hydrogen atoms at speeds up to — I30km/s 
identified by VM003) is not confirmed in the BJ007 anal- 
ysis. Velocities of — I30km/s are above the escape ve- 
locity of about ~42km/s (at 1 R p ). If observed, high 
velocity atoms must be escaping the planet. 

Following the same approach as above using the 
VM003 dataset, we evaluate the absorption seen during 
transit within the "blue" and the "red" sides of the wave- 
length domain defined by BJ007. In this domain, we find 
(9.8±1.8)% and (5.2±1.0)% absorption in the "blue" and 
the "red" sides, respectively. Therefore, the large spec- 
tral domain defined by BJ007 (including velocities up to 
±200 km/s) shows a significant absorption in both blue 
and red sides with, as found in VM003, an absorption 
stronger in the blue than in the red. 

Finally, although very little work has been published 
to explain these high velocities, they can be produced by 
radiation pressure from the intense Lyman-a flux of the 
nearby star (see the velocity diagram and the cometary 
shape of the escaping atoms in Fig. 3 of Vidal-Madjar & 
Lecavelier des Etangs (2004)). In which difference 
between blue and red absorptions is expected. 

4. THE LYMAN- a STELLAR VARIATIONS 

In previous works, corrections for stellar variations was 
done either using the wings of the line (VM003) or the 
observations before and after the transit (VM004) , while 
BJ007 correction was done by averaging the flux varia- 
tions by fitting a standard transit curve over the phase- 
folded data points (see Fig. 2 in BJ007). 

Contrary to the claim made by BJ007, the Lyman a 
variations of HD 209458 are not necessarily relatively 
large but normal for a quiet solar type star (Vidal- 
Madjar, 1975). In fact, when data are phase-folded, 
the apparent variations can be artificially enhanced due 
to the superimposition of observations made at different 
dates. Moreover, real stellar fluctuations are combined 
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TABLE 1 

Evaluated Lyman a absorption over various spectral domains 



Line core Intermediate Whole line 

Spectral domain 1215.15-1216.10A 1214.83-1216.43A 1210-1220A 

VM003 limits BJ007 limits 



Published absorption (%) 15 ± 4 8.9 ±2.1 5 ± 2 

(reference and data set) (VM003) (BJ007) (VM004) 

Absorptions using VM003 dataset a 15.1% ±4% 7.3% ± 2.0% 5.7% ±1.9% 

a Absorption evaluated using VM003 datasct 

and following the VM003 method with wings as flux reference. 
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Fig. 2.— Plot of the Solar and HD 209458 "Core"/ "Wing" ratios 
as a function of the total Lyman a flux ( "Core" + "Wing" ) . "Core" 
and "Wing" domains are defined in Fig. [T] Triangles are for Solar 
values. Circle and square symbols are for HD 209458 and corre- 
spond to the three HST visits A, B, and C, each visit including 
three successive HST orbits, namely 1, 2 and 3 (see VM003 for de- 
tails). The linear regression to the Solar values is shown along with 
the upper and lower limits (dotted lines). The HD 209458 values 
show basically two groups. In a first group, out-of-transit measure- 
ments Al, Bl and CI are within the limits of the Solar fluctuations. 
The second group consists of in-transit measurements (where A2, 
B3 and C3 are the closest mid-transit) that are well below the 
Solar fluctuations. The comparison between the dispersion within 
a group and the amplitude of the transit signature between the 
two groups, the absence of correlation between the HD 209458b 
"Core" / "Wing" ratios and the total Lyman a flux, both invali- 
date the suggestion that HD 209458 Lyman a activity corrupts the 
transit signature. 

with significant noise which is expected to be large be- 
cause of the limited number of photons in the short tem- 
poral bins. 

The Sun can be used as a proxy to study the possible 
variations of the flux and shape of the HD 209458 Ly- 
man a line. HD 209458 is a Solar type star (GOV) whose 
similarity with the Sun is demonstrated by the observa- 
tion of the chromospheric Call line profiles, which are 
relatively quiet for both stars. From observations ob- 
tained with SOHO (Solar and Heliocentric Observatory, 
Lemaire et al., 2002), we show that the core-to- wing ra- 
tio of the Solar Lyman a line varies by about ±8%, while 
it varies by about ±6% for a given total Lyman a flux 
(Fig. [5]). HD 209458 is expected to have the same be- 
havior and the 15% transit measurements is unlikely to 
be due to stellar variations. This is strengthened by the 
HD 209458 Lyman a measurements which present two 
different groups of core-to-wing flux ratios (Fig. [2]) : the 



ratios measured before the transit and the ratios during 
the planetary transit. Although these measurements arc 
obtained at various epochs, the dispersion within a given 
group is small and within the error bars, while the differ- 
ence between the measurements taken during the transit 
compared to the reference measurements taken before the 
transit is significant and larger than variations observed 
in the Sun. This behavior is unlikely to be coincidental 
and corresponds to the expected signal for a cloud of HI 
atoms in the environment of HD 209458b. 

Moreover, the dispersion of these individual measure- 
ments shows no correlation with the total Lyman a flux. 
This shows that the stellar variations are unlikely to be 
responsible for the observed Lyman a core-to-wings vari- 
ations. The core-to wing variations are thus more likely 
related to the planetary transit. 

Finally, it is extremely unlikely that stellar variations 
can mimic a transit light curve as seen in Fig. 2 of BJ007, 
when measurements arc phase-folded with the planetary 
orbital ephemerides. 

In summary, the stellar variations are taken into ac- 
count by VM003 by comparing the variations in the core 
to the variations in the wings of the line, and by BJ007 
by averaging a smooth transit light curve over a fluctuat- 
ing Lyman a stellar flux. The similarity of the resulting 
absorption depths, using the BJ007 and VM003 dataset, 
when evaluated over the same spectral domain (Table 1) 
shows that the Lyman a stellar variability does not cor- 
rupt the transit evaluation whatever the approach used 
for the data analysis. 

5. THE SIZE OF THE ABSORBING CLOUD AND THE 
ROCHE LOBE 

BJ007 compares the 8.9% absorption depth derived in 
his work to the absorption caused by an occulting disk 
with a size of ~ 4.08 Rp supposed to be the size of the 
Roche lobe as calculated using equations found in Gu 
et al. (2003). In the case of HD 209458b, a disk with a 
radius of ~ 4.08 Rp corresponds to an absorption depth 
of about 25% during transit. BJ007 thus concludes that 
the observed hydrogen atoms are inside the Roche lobe 
and cannot escape the planet. 

First, as discussed above, the BJ007 evaluation of 
~ 9% is only a fraction of the full HI absorption depth in 
the line core. Second, the formula of Gu et al. (2003) for 
computing the Roche lobe corresponds to the Lagrangian 
point LI between the star and the planet, i.e. the most 
distant Roche limit position relative to the planet (see 
Eq. B.8 and discussion in Lecavelier des Etangs, 2007). 
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However, the Roche lobe around the planet is not spher- 
ical but elongated toward the star (see e.g. Lecavelier 
des Etangs et al. 2004). In a transit configuration, the 
observed limit of the Roche lobe is in a direction perpen- 
dicular to the star-planet direction. In this perpendicular 
direction, the Roche lobe extension is about 2/3 of the 
extension toward the LI point. Therefore, it is more ap- 
propriate to use an average distance to the Roche lobe, 
which was given by VM003 to be ~ 2.7 Rp (Paczynski, 
1971). A filled Roche lobe corresponds to about 12% ab- 
sorption. This value is comparable to the HI observation 
of > 10% absorption depth in the line core. We can con- 
clude that HI atoms reach the Roche lobe or beyond, in 
agreement with the models of atmospheric escape. 

In addition, although observation of hydrogen atoms 
outside the Roche lobe is a direct evidence for escape, 
it is not a necessary condition. Even filling-up half a 
Roche lobe would imply escape rates large enough to 
significantly affect atmospheric structure (Lecavelier des 
Etangs et al. 2004). Thus, whether the hydrogen cloud 
actually fills up the Roche lobe or not, the large extension 
of the upper atmosphere (and possibly a cometary shape 
due to radiation pressure) shows that the atmosphere is 
indeed escaping. 

6. CONSEQUENCE: EVAPORATION IS CONFIRMED 

We have demonstrated that the BJ007 analysis is in 
agreement with HI escape from the HD 209458b atmo- 
sphere. Even if one considers that the absorption value 
given by BJ007 is correct, then either: i) this ~9% ab- 
sorption depth is taking place over the whole spectral 
range of ±200km/s, or ii) the ~9% absorption is the 
result of an unresolved larger absorption within a nar- 
rower wavelength range, produced by atoms below the 
escape velocity of about 42km/s. In the first case, HI 
atoms are detected to move at velocities much larger 
than the escape velocity. In the second case, the ab- 
sorption takes place over a narrow spectral range, ~5 
times narrower than the range considered in BJ007. In 
this narrow wavelength range, the absorption must be 
»12% and HI atoms must be present beyond the Roche 
lobe. In both alternatives, the result of the data analysis 
described in BJ007 shows that HI must be escaping the 
planet atmosphere. 

7. CONCLUSION 

BJ007 called into question the VM003 discovery of 
the atmospheric escape from the HD 209458b extrasolar 
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planet. BJ007 gives two main arguments : i) the absorp- 
tion depth is smaller than previously estimated and the 
Roche limit is not reached, and ii) the data analysis is 
corrupted by intrinsic stellar Lyman a variations. The 
first argument is not correct because the absorption is 
not taking place over the whole Lyman a line. The ab- 
sorption depth measurements depend on the considered 
wavelength range. In addition, the Roche lobe shape and 
orientation must be taken into account when comparing 
the size corresponding to these absorption depths with 
the size of the Roche lobe. The second objection is re- 
jected by recalling that VM003 corrected for the stellar 
variations using the overall stability of the line shape and 
the line wings as flux reference. 

Finally, despite the affirmation in the introduction of 
BJ007, works to transpose these observations into esti- 
mates of escape rate have been made early in VM003 
and recently in Schneiter et al. (2007). As stated by our 
anonymous referee, these "observational" estimates agree 
with the most recent and sophisticated models which 
all find rather good agreement with mass loss rates of 
~ 3 — 7 x 10 10 gs _1 . Therefore, at the orbital distance of 
HD 209458b, atmospheric escape does not strongly affect 
the evolution of the planets, in accord with the evolution- 
ary studies. 

Last, but not least, Ben Jaffel (2007) performed an 
independent, thorough and careful analysis of the best 
available data set for measuring the Lyman a transit ab- 
sorption depth. Despite a misinterpretation of the result- 
ing spectra, this work confirms that the signal detected 
in VM003 is indeed present and is not related to the 
data reduction process. Therefore, in the atmosphere of 
HD 209458b, HI appears to be the only species to have 
been detected with two independent datasets, and now 
with independent data analyses. The BJ007 data anal- 
ysis strengthens the detection of the extended escaping 
atmosphere of HD 209458b, in agreement with numerous 
theoretical studies of this phenomenon (e.g., Lammer et 
al. 2003; Lecavelier des Etangs et al. 2004; Yelle 2004, 
2006; Baraffe et al. 2005; Tian et al. 2005; Garcia-Munoz 
2006). 
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